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Abstract 

The influence of contacts on linear transport through a molecular wire attached to 
mesoscopic tubule leads is studied. It is shown that low dimensional leads, such as 
carbon nanotubes, in contrast to bulky electrodes, strongly affect transport prop- 
erties. By focusing on the specificity of the lead-wire contact, we show, in a fully 
analytical treatment, that the geometry of this hybrid system supports a mecha- 
nism of channel selection and a sum rule, which is a distinctive hallmark of the 
mesoscopic nature of the electrodes. 
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1 Introduction 

Future electronic miniaturization may enter a regime where devices are dom- 
inated by quantum mechanical laws and eventually reach the single-molecule 
scale [1] . Although molecular materials for electronics have already been real- 
ized [2], their implementation in real applications [3-5] still has to cope with 
challenges in utilization, synthesis, and assembly [6]. Concerning theoretical 
ideas and methods the problem is also two-sided: On the one hand many the- 
oretical ideas are already footed on past pioneering work, such as the proposal 
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of molecular rectification in 1974 [7] which was experimentally realized only 
20 years later [8]. On the other hand most of the conventional methods fre- 
quently employed for characterizing transport properties in microelectronic 
devices, such as the Boltzmann equation, can no longer be applied at the 
molecular scale. Here conductance properties have to be calculated by em- 
ploying full quantum mechanical approaches and by including the electronic 
structure of the molecules involved. 

Electron transport on the atomic and molecular scale became a topic of intense 
investigation since the invention of the scanning tunneling microscope (STM). 
More recently, studies of transmission properties of single molecular junctions 
contacted to metallic leads [9,10] have intensified the interest in the basic 
mechanisms of conduction across molecular bridges. The archetype of such a 
molecular device can be viewed as a donor and acceptor lead coupled by a 
molecule acting as a bridge. In such systems the traditional picture of electron 
transfer between donor and acceptor species is re-read in terms of a novel view 
in which a molecule can bear an electric current [11]. Molecular bridges have 
been realized out of single organic molecules [9,12], short DNA strands [10], 
but also as atomic wires [13-15]. Generally, contact effects alter the "intrinsic 
conductance of the molecules in such experiments and call for closer theoretical 
studies. 

In a parallel development the use of carbon nanotube (CNT) conductors has 
been the focus of intense experimental and theoretical activity as another 
promising direction for building blocks of molecular-scale circuits [16,17]. Car- 
bon nanotubes exhibit a wealth of properties depending on their diameter, 
orientation of graphene roll up, and on their topology, namely whether they 
consist of a single cylindrical surface (single-wall) or many surfaces (multi- 
wall) [18]. If carbon nanotubes are attached to other materials to build ele- 
ments of molecular circuits, the characterization of contacts [19,20] becomes 
again a fundamental issue. This problem arises also when a carbon nanotube 
is attached to another molecular wire, a single molecule or a molecular cluster 
with a privileged direction of the current flow. 

In the usual theoretical treatment of transport through molecular wires, the 
attached leads are approximated by a continuum of free or quasi-free states, 
mimicking the presence of large reservoirs provided by bulky electrodes. How- 
ever such an assumption may become inadequate when considering leads with 
lateral dimensions of the order of the bridged molecule, as for CNTs [19]. The 
latter have been recently used as wiring elements [17], as active devices [17,21], 
and, attached to scanning tunneling microscope (STM) tips, for enhancing 
their resolution [22,23]. With a similar arrangement the fine structure of a 
twinned DNA molecule has been observed [24]. However, CNT-STM images 
seem to strongly depend on the tip shape and nature of contact with the 
imaging substrate [25]. This calls for a better characterization of the contact 
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chemistry of such hybrid structures. 



This paper addresses the influence of the molecular wire-electrode contacts 
on the linear conductance when the spectral structure and the geometry of 
the electrodes plays an important role. This allows us to quantify to which 
extent mesoscopic leads may affect the conductance. Owing to the relevance 
of CNT-based devices, we focus on bridges between tubular leads. In previous 
density-functional-theory based treatments the conductance through systems 
such as a Cqq molecule in between two CNTs has been calculated with high 
accuracy [26]. These numerical approaches showed a strong sensitivity of the 
current on the system geometries and strength of the molecule-CNT couplings. 

In the present paper we focus on such contact effects. As a simple model for an 
atomic or molecular bridge we use a homogenous linear chain which enables 
us to derive analytical expressions for the conductance in a non-interacting 
electron approximation. In addition, we implicitly assume that no significant 
charge is transferred between the leads and the molecular bridge at equilibrium 
since this could lead to an electrostatic potential-induced inhomogeneity [27]. 
The latter may hold for an all-carbon [26] structure and makes it possible to 
investigate the properties of our model in the whole parameter space. The sys- 
tem exhibits distinct transport features depending on the number and strength 
of contacts between the molecular bridge and the interface as well as on the 
symmetry of the channel wave functions transverse to the transport direction. 
Our findings, which are common for leads with tube topology, are then stud- 
ied in detail for CNT leads (Fig. 1.) by analytically treating the single-particle 
Green function. In particular, we demonstrate on the one hand that config- 
urations with only one molecule-lead contact activated give rise to complex 
conductance spectra exhibiting quantum features of both the molecule and the 
electrodes; on the other hand multiple contacts provide a mechanism for trans- 
port channel selection leading to a regularization of the conductance, entirely 
provided by topological arguments. Channel selection particularly highlights 
the role of molecular resonant states by suppressing details assigned to the 
electrodes. 



2 System and Method 

In a tight-binding description, the hamiltonian of the entire system, H = 
-^tubcs H w [ ve -\- H coupling? reads 




a=L,R> w i re n a< n 'a 
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Fig. 1. Schematic representation of the molecular wire-carbon nanotube hybrid 
with single (bottom) and multiple (top) contacts, on-site energies e a=L > R ,wire are 
chosen to be zero. 

- E r mL (c^cf e + h. c .) - £ r mR (icf + h.c.) . 

m L <M c m R <A/ c 

Here, the matrix element „, = 5 n n i — 7/1 „, \ contains the on-site 
energy of each of the n w - 1TC = 1, . . . , N chain- atoms, e wnc , the orbital energy 
relative to that of the lead atoms, e L,R , and 7 L ' R , <-y wirc ; and T are nearest 
neighbour hopping terms between atoms of the left or right leads, molecular 
bridge, and the bridge/lead interface, respectively. Note that n Li R is a two- 
dimensional coordinate spanning the tube lattice. Summations over ?til and 
run over interfacial end-atoms of the leads. In general, there are M such 
atomic positions, defining the perimeter of the tube ends. The number of 
hybridization contacts between a tube and the bridge range between M c = 1 
(single contact case, SC) and M c = M (multi-contact case, MC). Typical real 
molecular wires are 7r-conjugate carbon chains with thyol end groups which 
in the present treatment are replaced by the linear chain model. 

Since the major results we present are not qualitatively affected by the use of 
more realistic quantum chemical models which take into account the precise 
structure and properties of the molecular bridge and of the attached leads, we 
keep the description of our problem at the level of the tight-binding model. In 
order to highlight the topological properties of tubular leads, we first study 
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the simplest case in which periodic boundary conditions are imposed on semi- 
infinite square lattice stripes, with the cuts parallel to the lattice bonds. We 
call this electrode specie square lattice tubes (SLT). In the case of CNT, when 
the graphene honeycomb lattice is rolled along the lattice bonds such that I 
hexagons are transversally wrapped, an armchair single wall (£, £) nanotube 
is obtained, and M = 21. According to Eq. (1) CNT are then described at 
the single-band tight-binding level for tt orbitals. This is equivalent to as- 
sume that among the four valence carbon orbitals no interaction between the 
a (2s and 2p x y ) and n orbital is significant because of their different symme- 
tries. The fourth electron, a p z orbital, determines the electronic properties 
which can be calculated by means of a tight-binding treatment, on the same 
level as we treat the molecular bridge. There are two such electrons per unit 
cell in a honeycomb structure, the 7r and ir* band, rendering the electronic 
properties of the material interesting, i.e. it can be a priori either metallic or 
semiconducting [28]. 

We study quantum transport in the framework of the Landauer theory [29] 
which relates the conductance of the system in the linear response regime to 
an independent-electron scattering problem [30]. The electron wavefunction 
is assumed to extend coherently across the whole device. The two-terminal 
conductance g at zero temperature is simply proportional to the total trans- 
mittance, T(E-p), for injected electrons at the Fermi energy Ep: 

g = (2e 2 /h) T(E F ) . (2) 

The factor two accounts for spin degeneracy. The transmission function can 
be calculated from the knowledge of the molecular energy levels, the nature 
and the geometry of the contacts. One can see this by expressing the Green 
function matrix of the full problem, G = G wirc +S L + S R , in terms of the 
bare wire Green function and the self-energy correction due to the presence of 
the leads. Making use of the Fisher-Lee relation [31] one can finally write 

T(E) = 4 Tr {A l (E)G(E)A r (E)&(E)} , (3) 
where 

A a (E) = ~ (E"(s) - E at (*)) 

For the system under investigation where only the first and last atom of the 
chain is coupled to the leads, the formula for the transmission simplifies to 

T(E) = AA L (E)A R (E)\G 1N (E)\ 2 , (5) 



(4) 
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where the spectral densities A L and A R are the only non-zero elements ^A L J 

and (A R ) , respectively, of the matrices A. The matrix element A L ( R ) is 
the spectral density of the left (right) lead. It is related to the semi-infinite 
lead Green function matrix Cj L ( R \ It is minus the imaginary part of the lead 
self-energies (per spin), 



s* = A*-iA Q = £ r m r* m , a g« am , a , (6) 

m a ,m' a 

with a = L, R. Owing to the causality of self-energy, its real part A can be 
entirely derived from the knowledge of A via a Hilbert transform. 

The rhs of Eq. (5) coincides with formulas used to describe electron transfer 
in molecular systems [32] . The above relationship between the Landauer scat- 
tering matrix formalism on the one side and transfer hamiltonian approaches 
on the other side has been worked out in the recent past [11,33] showing de 
facto their equivalence. This enables us to make use of the formulas from a 
Bardeen-type picture in terms of spectral densities, which is often convenient 
for an understanding and analysis of results obtained. 



3 Molecular Green function 



One has to calculate from the N x N matrix, G wirc 1 = E + i0 + — H wiic , the 
Green function matrix element Gin needed in Eq. (5). This matrix element 
refers to the two ends of the A^-atom-molecule. Its computation requires an 
N x N matrix inversion. Since only the molecular-end on-site energies are 
perturbed by the interaction with the leads via the self-energies S Q , some 
general conclusions can already be drawn without an explicit computation of 
Gin, namely one can write [34] 



G 



wire 
IN 



Gin — y — ■ . ^ T — . — . To - (7) 



£ L G7r) (1 - T, r G n 1 n) - S L E R (G 



wire 
IN 



The interaction with the leads dresses, via the self-energy E a , the bare molec- 
ular wire Green function element GY§°. The latter can be calculated analyt- 
ically in the case of a homogeneous wire {e^ ire = e wire , 7™ e n / \ = 7 wirc ). In 
fact, upon projecting on the N dimensional molecular wire basis, the deter- 
minant of the bare molecular Green matrix factorizes into a dimensionless 
function of only the number of chain atoms, and of the ratio rj wac = (E — 
e wim ) I (27 wirc ). This leads to a closed form for the molecular contribution to the 
conductance. Namely, one can easily check that GJ§ C = <y wircN 1 det (G v 



"fwire ' 
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(7 wirc ) 1 £o/£n, and G™ re = Gffi = 7 wirc ^n-i/Zn, where the exact form of 
£ reads: 



/ , \ N+l / , \ N+l 

(v wirc ) = [v wivc + vV ire ) - ij - (v wirc - vV irc ) - 1) • 

After some algebra one finds that £ possesses the recursive property £,n£,n-2 
£jv-i ~~ £o> which leads us to re-write Eq. (7) as 



CJV ~~ : 1 : — CJV-1 H : — 9"SA r -2- 

I . wire. - v wire / ' 1 _,wire z 



7 wire GlAr V 7 wirc Y 



In other words, the inverse Green function matrix element connecting left and 
right leads can be written as a sum of terms, representing the inverse of the 
bare Green function matrix elements for a wire of N, N — 1, and N — 2 atoms. 
In the limit of weak contact coupling, the behavior of the Gin element is 
dominated by £at, leading to iV transmission resonances in the conductance 
of unit height. Nevertheless, if the effective coupling between the molecule 
and the lead is much larger than 7 wire , £ N _ 2 will become the dominant term. 
As a consequence the conductance spectrum is effectively that of an (N — 2)- 
atomic wire [35]. The resonant behavior inside the wire band (|?7 wirc | < 1) 
and its modification due to the lead coupling is easily understood by writing 
the transmission as T = 45 2 sin 2 (0)/Z> where the denominator T> takes the 
following compact exact form valid for all N > 1: 



V = (sin(A + 1)0 - (<5 2 - A 2 ) sin(A - 1)0 - 2A sin A0 
+ 4 (5 sin A0 - A<5sin(A - 1)0) 2 . 

Here a = A — i5 = S/7 wire is the self-energy of the leads (for simplicity assumed 
to be equal) normalized by the wire hopping. The parameter 



-wire _ /„ w ire 



= cos" 1 7T re = \ In V , = (9) 

£ rewire _|_ ,1 j-jwire- 1 J 

is real in the wire band giving rise to resonances for injected electrons match- 
ing the wire eigenenergies. Outside the wire band is pure imaginary (sin 
functions are effectively sinh functions), and the transmission has a power law 
dependence on energy and an exponential dependence on the wire length, that 
is 

T ~ |2r/ wirc |- 2Ar for |^ wiro | > 1, (10) 



,2 
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in agreement with previous results [36]. This analytic expression for the trans- 
mission provides the generalization of existing results [32,33,37] to the case 
with non-vanishing real part of the self-energies. The density of states M = 
— Im Tr {G} /ir can also be written in a closed analytical form. One can there- 
fore take advantage of the fact that, due to the wire homogeneity, all the 
diagonal elements except the first and the last one coincide, 



Kjrhh — 

1<k<N 7 wire 6v - 2(7^-i + a^ 2 



By using the parametrization (9) one can easily recast the density of states 
into the compact form 



M = x 

N sin m - 2(N - l)a sin(iV - 1)0 + (N - 2)a 2 sin(iV - 2)0 



Im 



sin(iV + 1)0 - 2a sin N$ + a 2 sin(iV - 1)0 



4 Electrode self-energy 



In calculating the spectral function, we make use of the assumption of iden- 
tical left and right leads and drop the self-energy indices in Eq. (6). Since 
the Hamiltonian is discrete, we can write the lattice Green function G = 
(E + i0 + — H)" 1 in matrix form by rearranging the two dimensional n lattice 
coordinate in Eq. (1). We assume the x direction to be parallel to the tubes 
(and to the transport direction) and y to be the finite transverse coordinate 
(see Fig. 1). The latter is curvilinear with n y spanning M sites with periodic 
boundary conditions. 

The lattice representation of the lead Green function is needed in the calcula- 
tion of the self-energy contribution. It can generally be written by projecting 
the Green operator onto the localized state basis, i^k x ,k y {n x = border, n y ) = 
Xk x <Pk y (n y ), of the semi-infinite lead: 



Gn v n> (E) 



E 



E + i0 + - H 



-i 



n, 



Xk x (pk y (n y )xljt y (n' y 



E + i0+ - E l 



(11) 



kx iky 
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4-1 One- dimensional electrodes: Newns model 



We first recall the particular case of linear chain electrodes (onto which we 
will map our system due to the validity of a channel selection). For such 
a simple model the dispersion relation as a function of the lattice on-site 
energies e, of the hopping terms 7, and of the lattice spacing a is simply given 
by E = e — 27 cos k x a. As a result the surface Green function for the semi- 
infinite chain is obtained by inserting the wave function at the lead origin 



Xk x — \J2/Tism.k x a in the defining expression (11) and transforming the only 
sum over momenta into an integral, due to the infinite system size. Thus, 



sin 2 k r a e ^ x {E)a 



-it /a 



(Ik 

x E + i0 + — e + 27 cos k x a —7 



where we solved the integral according to Refs. [38] and we made use of the 
dispersion relation. The resulting spectral density, given by Eq. (6), is the 
semi-elliptical local density of states (LDOS) as obtained by Newns in his 
theory of chemisorption [39]: 



r 

Newns ' 



2 



fo) = ^>/w e(i-H). 

r e fj is the strength of the single contact between the molecule and the semi- 
infinite one-dimensional leads, 77 = (E—e) / (27) is the band-normalized energy, 
and G the Heaviside function. The real part of the self-energy, responsible for 
shifting the molecular resonances, is simply proportional to cos k x a and thus 
linear in energy. Its full dependence on energy is given by Hilbert transforming 
A, and it reads 



r 2 

7 



(v - vV-i (e (v - 1) - g (-rj - 1))) . 



4-2 Square lattice tubular electrodes 



Square lattice leads are characterized by periodic boundary conditions per- 
pendicularly to the lead direction. Transverse momentum quantization leads 
to k J y a = 2ft j/M (with < j < M). The surface Green function for such a 
system can be written as 
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ir/a 

Sill 



^ V* v ; nM J * E + iQ+ _ g + 2 cog fc + 2 7 cos a 

fc « -it /a 

1 M_1 

= T7 M&(E)^(n' y ), (12) 



M . 

j=0 



where & (E) = —e lk *( E ) a / ^ has been obtained by solving an integral formally 
equivalent to the linear chain case and using the dispersion relation 



E = e-2-y (cos k y a +cos k 3 x (E)a) . (13) 

The transverse profile of the wave function is given by ipj(n y ) = exp(ik 3 y n y a) . 
Note that the wave function is obtained by a further normalization, namely 
= if /(Ma) 1 / 2 . 

The self-energy finally reads as a sum of weighted longitudinal wave function 
profiles 

i M-l 

e = m £ & (E) Vj/M [r] , 

1V1 j=0 
where the weight 



Vj/M [r] 



M-l 

E r m ^(m; 

m=0 



(14) 



is the contact-averaged transverse wave function. Depending on the contact 
geometry one has to specify the distribution of the T m contacts to calculate 
the weight rj and thus the self-energy. Note that rjr.) is formally the square 
modulus of the Fourier series of rV.); thus the zero-mode transverse momentum 
state j — contributes to r\ with the square of the mean contact strength. 
Due to the geometry of the lead surface, it is reasonable to assume a uniform 
distribution of contacts between the molecular wire and the electrodes. For 



contacts of equal strength T r 



active on M c < M sites, we obtain 



a modulation for the contributing channels governed by 



fr M \ r2 m sine 2 (irjMjM) 

vj/M r ef r, m c ) = r cff M c , 2 , 

sine (itj/M) 

where sinc(x 7^ 0) = sin(x)/x, and sinc(x = 0) = 1. One can decompose 
the spectral density into a sum over the spectral densities of each state j. 
Namely A = A (0) Y,j Wj(E), with A<°> = T^MJ^M). The channel weights 
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Fig. 2. Spectral weights Wj(E) plotted for different contact values M c ; a nonlinear 
grey level scale is used with black corresponding to and white to 1; small weights 
are amplified for better visualization. In every panel the horizontal axis represents 
the normalized energy —2 < {E — e)/{2^) < 2, and the vertical axis the normalized 
wave number < j/M < 1. Note that in the mesoscopic limit, M c /M < 1, the 
states j i / match the nodes of Aj(E): only the zero-transverse momentum j = 
contributes to transport (a better resolution figure is available upon request). 

are obviously independent upon rotation of the interfacial coupling position 
as S itself is. 

In Fig. 2., the weights Wj(E) are visualized for different contact values 1 < 
M c < 6. 

For the case M = M c the contributions from all states are suppressed except 
the state with zero transverse momentum, which is the outcome of the sum 
rule (14). That is, rj = F^gMSj^- Thus the configuration with all contacts of 
the tube ends coupled to the molecule with strength T c g/\^M is equivalent 
to the case of a single contact with strength T c s to a one-dimensional lead 
(previous section). Moreover a scaling law is found for E, and a fortiori for 
the conductance given by g = g [T ^MA , where f is the local contact strength. 

In Fig. 3, A is displayed as a function of energy, lead diameters and active 
contacts. As easily visible, it is only for values M c of the order of the avail- 
able contacts M that the mesoscopic nature of the scattering channels enter 
the spectral density. The larger the tube diameter the lower is the number of 
contacts which are needed to reach a MC-like spectral density. This observa- 
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Fig. 3. The normalized spectral density A7/r^ ff plotted for different maxi- 
mum values M c as in Fig. 2.; the x-axes represent the normalized energies 
—2 < (E — e)/(2-y) < 2, the y-axes label the dimensionality of the leads (num- 
ber of possible contacts) M c < M < 20 (a better resolution figure is available upon 
request). 

tion justifies the use of the one- dimensional Newns model for leads of lateral 
dimension much larger than the contacted molecule but also shows the limit 
of this approach when dealing with quasi-onedimensional leads. It remains to 
investigate to which extent the results obtained so far can be generalized to 
realistic quasi-onedimensional structures such as CNT. 



4-3 Carbon nanotube electrodes 



When the armchair (£, £) CNT topology is imposed the number of carbon sites 
at the interface is M = 2£. The eigenvalues of the tight-binding hamiltonian 
(1), 



E± (kl, j) = e ± 7^ 1 + 4 cos cos yY~J + 4 cos 




(15) 



are obtained in a basis set given by symmetric (+) and antisymmetric (— ) 
site configurations of the graphene bipartite lattice, corresponding to ir and 
7r* orbitals respectively [28,40]. The longitudinal momentum is restricted to 
the Brillouin zone, —it < k J x a < it, and the transverse wave number 1 < j < 2£ 
labels 4£ bands, as many as the number of atoms in the unit cell of a (£, £) 
CNT. The two bands corresponding to j = £ are singly degenerate. They are 
responsible for the metallic character of armchair carbon nanotubes (these two 
bands cross at the Fermi level E = e for k^a = ±2tt/3). Also the two outermost 
bands corresponding to j = 2£ are singly degenerate while the other remaining 
(4£ — 4) bands are collected in (2£ — 2) doubly-degenerate dispersion curves. 
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The single-particle Green function in a lattice representation for two sites 
belonging to the same sub-lattice can be still written as in Eq. (12) as 



1 2e 

$n y ,n> (E) = - £ <P, W & (E) $ (n' y ) , (16) 

where <fij(n y ) = exp(ik y n y a) , with k y a = nj/£, and 1 < j < 2£. Note that in 
Eq. (16), n y and n' y should be either even or odd (that is they should belong 
to the same sublattice). The semi- infinite longitudinal Green function is given 
by 

&(E) = — T T dk j sin2 ( fc >/ 2 ) 

The integral can be worked out analytically by extending k x to the complex 
plane and adding cross-cancelling paths (parallel to the imaginary axis) along 
the semi-infinite rectangle in the half plane lmk x > and based on the interval 
between —it /a and ir/a. The closing path parallel to the real axis gives a 
real contribution linear in energy. This generalizes the approach by Ferreira 
et al. [41], recently adopted for obtaining an analytical expression for the 
diagonal Green function of infinite achiral tubes, to the case of semi-infinite 
CNTs. The determination of the poles inside the integration contour, given 
by 



cos 



1 

-- cos 
2 



2\\, 2 7 



sin 



J7T 



allows for the calculation of the residues and thus of the surface Green function. 
One finds 



& (E) 



1 E- 

27 27 



sin 



1 + i- 



sin 



J7T 



(17) 



where the choice of the contributing pole through the branch parameter /3* = 
sign (E — e) has to be taken into account. The LDOS, obtained from the 
imaginary part of the surface Green function after Eq. (17) is plugged into 
Eq. (16), is shown in Fig. 4. It clearly differs from the LDOS of an infinite 
CNT as depicted for comparison in the right panel. As for the case of the 
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Fig. 4. Left panel: the normalized spectral density for a semi-infinite (£,£) CNT 
lead in the SC configuration; it corresponds to the LDOS at any atom site at the 
cut of the CNT lead. For comparison the dispersion relation and the LDOS of an 
infinite (£, £) CNT are shown in the middle and right panel respectively. Solid lines in 
the dispersion relation panel indicate doubly degenerate bands, dashed lines singly 
degenerate bands. Here £ = 10, and on-site energies and hopping terms refer to 
a = L, R-leads. 



SLT the pinning of the longitudinal wave function at the surface of the semi- 
infinite systems cancels all border zone anomalies when q±a matches multiples 
of 2tt. In infinite SLTs these states are the only resonant states (van Hove 
singularities) so that the surface LDOS of a semi-infinite SLT never diverges 
(as it is shown in the left panel of Fig. 5). On the contrary, in CNTs there are 
states with zero group velocity outside the border zone which are responsible 
for the singularities of the spectral density of semi-infinite CNTs (left panel 
of Fig. 4). 

The self-energy for a CNT lead is more complicate than the one for a SLT 
owing to the missing equivalence of the sites belonging to the two different sub- 
lattices. However, since the longitudinal part of the Green function, Eq. (17), is 
the same for all diagonal and off-diagonal terms of the surface Green function, 
the self-energy can still be cast into the form 
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SLT, M — 10 (5,5) CNT, M — 10 




Fig. 5. The normalized spectral density as a function of energy and active contacts is 
plotted for M = 10 possible atomic contacts available; on-site energies and hopping 
terms refer to a = L, R- leads. The right panel illustrate numerical results after 
Ref. [42] in full agreement with the analytics showed in the text. 



1 21 ~. 
However, for the calculation of 



it 

Vj/ifi] = T m<Pj(m) 



one has to specify the sub-lattice components of the transverse wave func- 
tion and whether they belong to a bonding or anti-bonding molecular state. 
Again the distribution of the T m contacts is needed in oder to calculate the 
weight r\ and thus the self-energy. Eq. (18) simplifies considerably in the SC 
case: r\ = T 2 . Since t] is uniform in j the self-energy is simply proportional to 
the diagonal semi-infinite Green function and, as a consequence, the spectral 
density is proportional to the local density of states (Fig. 4.). The MC case 
(r m = r e fj / \p2f) is also easily tractable leading to a sum rule over the possible 
conducting channels. However a direct proof is provided by the intuitive con- 
sideration that only the 7r-bonding state can contribute to the MC spectral 
density (all the other states have a non-constant spatial modulation provided 
e.g. in Ref. [43]). Following our notation the 7r-bonding state corresponds to 
j = £. Fig. (5) shows the spectral density in the intermediate regime between 
the SC and MC limits. The two different lead lattice structures carry the same 
physical information only in the MC limit case. 
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5 Discussion and concluding remarks 

It is interesting to recollect the results for the MC-spectral density, Amc ; hi all 
three lead models considered so far. Using the dispersion relation (13) and the 
surface Green function in Eq. (12) the spectral density for SLT leads coincides 
formally with the Newns LDOS with an energy shift (r/) = A Ncwns (rj + 1). 
For CNTs Eqs. (15-17) lead to Ag^ T (r/) = A Newns (77 + §). From the above 
discussion it is clear that the multiple contact configuration suppresses features 
associated with the two-dimensional character of tubular leads, apart from an 
energy shift. In our model, the latter is the only remnant the system preserves 
from the transverse momentum component. 

In contrast, the SC case is strongly dependent on the lead underlying structure. 
The spectral density for a single contact, Age, reduces to the LDOS in the 
lead at the point where the molecular wire is contacted with strength T, 

A sc = vrr 2 LDOS. 

It is, in particular, in the SC scenario that the conductance of the molecular 
wire might be strongly affected by the nature of the leads [35] . Nevertheless, 
once the nature of the contact can be inferred, one can think to cure the 
spurious insertion in the conductance by filtering out the contribution of the 
leads from the molecular resonances. For instance, in the CNT-enhanced STM 
tips [22] the improved resolution images can be cleaned by de-convolving them 
using model assumptions for the leads and their contact geometry. 

Another significant consequence of the peculiar contact dependence of the 
spectral density is the possibility to understand the influence of the mesoscopic 
character of the leads. In the limit of large M (at fixed M c ), rj looses its 
granularity being sampled by many more states compared to its nodes, whereas 
for M c /M < 1 an increasing number of nodes matches the decreasing number 
of states. This determines a reduction in the self-energy, and thus in the width 
of the molecular resonances, highlighting the quantum features of the wire. 
The latter result in a quite striking behavior for CNTs because of the band 
anomalies outside of the border zone which strongly determine the resonant 
behavior of the spectral density. 

To conclude, we have shown that novel features are expected to arise in the 
conductance of a molecular wire connected to nanotube leads. The commonly 
used approximation of a pure imaginary, flat, wide band self-energy is not 
valid when employing tubular leads. Nevertheless, the conductance of a homo- 
geneous molecular wire still possesses an analytical form in the entire regime 
of the wire parameters and allows for the insertions of a nonvanishing real 
self-energy, necessarily arising when considering nanotube leads. By tailoring 
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the geometry and dimensionality of the contacts, it is possible to perform 
a channel selection. In the MC limit the conductance becomes independent 
of the lattice structure of the tubular electrodes, transport is dominated by 
topology properties and is effectively one-dimensional. Furthermore, the con- 
ductance obeys a universal scaling law in the multiple contact configuration. 
We would further like to stress that the derived analytical expression for the 
semi-infinite CNT self-energy allows for a full analytical treatment of the lin- 
ear conductance problem. The possibility to handle an exact expression of 
the semi-infinite CNT Green function may serve as a first step in analytical 
treatments of more complex carbon based molecular structures such as T- or 
Y-junctions [44]. 
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